rm(list=ls())
library(pwr)

# power analysis (p1 is observed response rate for Turkish constituents)
power.prop.test(n=1522/2, p1=0.6325459, p2=NULL, sig.level = .05, power=.8) # .8 power would be .07 effect
